2-Dimensional Polymers Confined in a Strip 
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Single two dimensional polymers confined to a strip are studied by Monte Carlo simulations. They 
are described by A'^-step self-avoiding random walks on a square lattice between two parallel hard 
walls with distance 1 <C -D <S (i/ = 3/4 is the Flory exponent). For the simulations we employ 
the pruned-enriched-Rosenbluth method (PERM) with Markovian anticipation. We measure the 
densities of monomers and of end points as functions of the distance from the walls, the longitudinal 
extent of the chain, and the forces exerted on the walls. Their scaling with D and the universal 
ratio between force and monomer density at the wall are compared to theoretical predictions. 



I. INTRODUCTION 

The behaviour of flexible polymers in a good solvent 
confined to different geometries and in the presence of 
walls or other obstacles have been studied for many 
years 0, . A particularly simple geometry is the space 
between two parallel walls. For simplicity we shall only 
discuss here the case of walls without energetic effects, 
i.e. the walls play a purely geometric role. 

An important theoretical prediction is that near such a 
wall the monomer density profile increases as z^^'^ , where 
z is the distance from the wall {z ^ D, and D is the 
width between the two parallel walls) and is the Flory 
exponent JJ . This should hold in any dimension of space 
d. On the other hand it is intuitively obvious that the 
force exerted by the polymer onto the wall is propor- 
tional to the monomer density near the wall. The ratio 
between the two can be expressed in terms of a univer- 
sal amplitude ratio which has been calculated by Eisen- 
riegler^3 (using conformal invariance results of Cardy 
et al. Pi) in d = 2, and in any d < 4 by Eisenriegler Q 
by means of an e-expansion. Attempts to verify these 
detailed predictions by Monte Carlo simulations in three 
dimensions |3i H have so far been without very con- 
vincing results. As far as we know, no attempt was made 
yet to verify them in d = 2, and that is where the present 
paper sets in. 

We study single polymer chains confined to a 2-d 
strip. They are described by self-avoiding random walks 
(SAWs) of TV steps on a square lattice between two hard 
walls with distance D as shown in Fig. ^ More pre- 
cisely, monomers are supposed to sit on lattice sites and 
D is the number of rows accessible to monomers, i.e. the 
walls are placed a.t y — and at y — D + 1, and the 
monomers can be at y = 1, . . . D. We only consider the 
case where the Flory radius of a free chain of length N, 
Rp K N" , is much larger than D. When using a chain 
growth algorithm, the polymer has then to grow , after a 
short initial phase of ~ D^^'^ steps, in either the positive 
or negative x-direction without possibility to change its 
orientation. This allows us to use an additional wall at 
X < xo which forces all chains to grow into the positive 
x-direction. For N '3> D^^^ this will essentially reduce 
the partition sum by a constant factor, without affecting 



any of the scaling laws or any of the detailed compar- 
isons with theoretical predictions. On the other hand, it 
simplifies the subsequent discussion. 

The force exerted onto the wall is most straightfor- 
wardly expressed in terms of the work done when moving 
one of the walls, i.e. by the dependence of the free energy 
- and thus also of the partition sum ~ on D, 
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where we have introduced a dummy temperature T which 
can take any positive value. The partition sum Zjy is just 
the number of TV-step SAWs in the strip starting from a 
given x, but summed over all values 1 < ya D oi the 
^/-component of the starting point. 

The partition sum of a free SAW in infinite volume 
scales for iV — > oo as Zm = ji'^ N^~^ with /ioo being the 
critical fugacity per monomer, and with 7 = 43/32 being 
a universal exponent. In contrast, the partition function 
on a strip scales as 
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without the power correction and with scaling for 
large D as 
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with a being a non-universal amplitude. The force per 
monomer is then obtained as 



/ = F/N - keT D 
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Standard fixed-length Markov Chain Monte Carlo sim- 
ulations do not give estimates of the partition sum or of 
the free energy, so that Eqs. H1I4I) cannot be used di- 
rectly. This has led to algorithms specifically designed 
for estimation of forces [lCll | , but emplo ying the pruned- 



enriched-Rosenbluth method (PERM) one can use 
Eqs. (|1I4|I directly. 

Using PERM with fc-step Markovian anticipation 0, 
IT^ IT^ , we measured the partition sum directly and es- 
timated the dependence of the monomer fugacity on the 
width D. In the same simulations also the monomer den- 
sity profile, the end-to-end distance along the strip, and 
the density profile of chain ends are measured. 
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FIG. 1: Schematic drawing of a polymer chain growing in- 
side a strip, and with an additional wall added at x = 0. 
Monomers are only allowed at lattice sites a; > and 1 < y < 
D. 



Details of the simulation method are given in the next 
section, results and their comparison with theoretical 
predictions are discussed in Sec. III. 



II. ALGORITHM: PERM WITH fc-STEP 
MARKOVIAN ANTICIPATION 

PERM is a chain growth algorithm with popula- 
tion control. Polymer chains are built like random walks 
by adding one monomer at each step. We use a Rosen- 
bluth like bias for self-avoidance. As usual, this bias 
is compensated by a weight factor 0|, i.e. each sam- 
ple configuration should be given a weight. But actu- 
ally we use a stronger bias which in addition suppresses 
dense configurations and samples more finely relatively 
open chain configurations, called Markovian anticipation 
in 0. 

In fc-step Markovian anticipation, we let the addi- 
tional bias in the next step depend on the last k steps 
made before. Let us denote the 2c? directions on a d- 
dimensional hypercubic lattice by s = 0, . . . ,2d — 1. 
All possible (fc + l)-step configurations {k previous steps 
i = —1, —2, . . . — k plus one future step i = 0) are then 
indexed by 
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In addition, choose an integer m « 100 (the precise value 
is not important). Either during an auxiliary run or dur- 
ing the early stages of the present run, we obtained a his- 
togram Hm such that iJm(S) is the accumulated weight 
at chain length n -f to of those chains which had con- 
figuration S during steps n — k,n — k + 1, . . .n. Since 
we sample uniformly, Hq{S) is independent of S. Thus 
iJm(S) with TO > indicates how "successful" is a config- 
uration S after to more steps. In importance sampling we 
want each chosen direction to have in average the same 
later success. Therefore we choose the next step with 
probability 



p(so|s) 



Hmis, So) 



(6) 



In our simulations we choose fc = 9 and to = 100. We 
accumulate contributions to H only for n + m > 300, and 



we apply Eq. © only for chain lengths > k (for chain 
lengths < k there is not yet enough history to condition 
upon) . 

Accumulating the histogram only for n ^ 1 is sug- 
gested by the fact that only for large n the anisotropic 
bias is fully developed This anisotropy makes the 

histogram strongly dependent on D. We found that us- 
ing for all D only the histogram obtained for free chains, 
i.e. for D = oo, gives nearly the same efficiency. This is 
quite different from the case where the anisotropy is not 
due to geometry, but is due to stretching of the polymer. 
In the latter case, simulations with markovian anticipa- 
tion become much more efficient with increased stretch- 
ing This is not the case for the present problem 
where the anisotropy is due to geometric constraints, 
for which PERM with Markovian anticipation is how- 
ever still the most efficient known simulation method by 
far. 



III. RESULTS 

Before presenting our results, let us stress that we have 
several possibilities for checking our algorithm. For very 
large D we can compare our estimates ofu with the very 
precise estimate ^oo = 0.37905228 [Ullll- For D < 11 
we can compare with exact transfer matrix calculations 
of 19]. And for D < 2 we can even solve the problem 
analytically. 

For D = 1, the polymer can only grow in a straight 
configuration, giving fii^D=i) — 1- For D — 2, each step 
can be either up (u), down (d), or to the right (r). After 
an 'u' or 'd' move, the next step has to be 'r', while any 
move is possible after 'r'. Lumping moves 'u' and 'd' 
together into a vertical move ('v'), we see that the set of 
all possible configurations forms a regular language 
with the associated graph shown in Fig. |21 The partition 
sum for chains of length TV is just twice the {N + l)-st 
Fibonacci number. 



ZnID = 2) - 2Fn+i (TV > 1) , 



(7) 



where Fq = Fi — 1 and F/v = ^V-i + Fn^2- Thus the 
critical monomer fugacity for D = 2 is the inverse of the 
golden mean. 



MD=2 = g 



= 0.61803. 



(8) 



In this case, we can also show that markovian antici- 
pation gives the optimal bias. Markovian anticipation 
corresponds in this case to : = 5 : 1 if a vertical 
move is allowed, and pr : Pv — ^ '■ else. This choice 
leads to weights which oscillate between two values, thus 
no population control (pruning/cloning) is needed f2]| . 
All this is verified in our simulations, which serves thus 
as a test for our algorithms. 

For D > 2, we used simulations. We simulated strip 
widths up to D = 320 and chain length between N = 
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FIG. 2: Graph accepting the regular grammar of vertical ('v') 
and right ('r') moves for D = 2. The node labelled 'A' is 
the start node |20|| . If one wants to distinguish also between 
'u' and 'd' moves, the graph is somewhat more complicated 
and contains also a transient part. This is skipped here for 
simplicity. 
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FIG. 3; Log- log plot of /xu — /loo against D. The dashed line is 
MD -Mcx) = 0.737/)"^/'^ with ly = 3/4, as predicted by Eq. @ 



3000 (for £) = 2) and iV = 125, 000 (for D = 320). Crit- 
ical fugacities are determined by plotting Z^^^ against 
log N and demanding that these curves become horizon- 
tal for large N. Results are shown in Fig. |3| where we 
plot " Moo with ^oo = 0.37905228 18]. They are 
in perfect agreement with the theoretical prediction of 
Eq. lO, and provide the estimate a = 0.7365 ± 0.0007. 
In addition, fijj can be compared for _D = 3 to 12 with 
the transfer matrix results of For all D, the values 
agree for at least six digits. 

As we said in the introduction, we expect the force / 
onto the wall to be proportional to the monomer den- 
sity p{y) near the wall. The precise relation is given by 
Eisenriegler Q 
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Here the non-universal amplitude k relates the end-to- 
end distance of a free SAW to the chain length, k — 
rI'" /N = 0.5297± 0.0002 jH for the square lattice. On 
the other hand, _B is a universal number. For ideal chains 
B = 2, while for chains with excluded volume in 4 — e 
dimensions one has B « 2(1 - bie) with h = 0.075 Q. 
The latter is of course of dubious value for d = 2 (where 
it would give B « 1.7), but conformal invariance leads to 
the supposedly exact value B = 2.01 in d = 2 HiQ. 
The monomer density is normalized such that 
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FIG. 4: (a) Rescaled values of the monomer density, {D + 
1) p{y) against ^ — y/{D + 1). Also plotted is the function 
/p(C) = 10.38 iai-Or^'- 

(b) The same values (for D < 160), but divided by fp{^)- In 
this panel we do not display our data for the largest lattices, 
since they are too noisy and would just blur the picture. They 
do however show the same trend as the data for D < 160. 

(c) The data for D = 128 plotted against a modified scaling 
variable, ~ {y — S)/{D + 1 — 25), divided by fpi^s), for 
three different values of S. 



Y^y=iP{y) ~ 1- According to Eq. ^ it should scale 
as y^/'^ near the walls (this holds in any dimension) with 
u = 3/4. Surprisingly, we found that the simplest ansatz 
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generalizing this power law to all y G [1, -D] 
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FIG. 5: End-to-end distance divided by A'', {x)/N, plotted 
against A*' for various values of D. 
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FIG. 6: log-log plots of Aa; and of the boundary density ph 
versus D. The dashed line is 0.915 D~^'^ and the solid line 
is 10.75 D-^. 
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FIG. 7; Rescaled values of the probability pendiy) that the 
chain end is at the distance y from a wall against ^ = y/{D + 
1). The topmost (interrupted) line is the function gp{^) = 
2.85 (e(l-C))''''''*'- 
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with A = r(14/3)/r2(7/3) = 10.38 gives already an ex- 
cellent fit (Fig. 2^^) in the limit D ^ oo. There are small 
deviations invisible in Fig. ^ but clearly seen when plot- 
ting p{y)/fp{C) (see Fig. but they seem to vanish 
slowly in the limit D ^ oo (see Fig. 0]d) . 

Assuming the latter, i.e. assuming that Eq. ifTUI) 
becomes exact asymptotically, the universal number B 
could be estimated by using Eqs. iQJ, Q, ®, and © 
which together would give 
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Inserting the above numbers gives B — 2.122 ± 0.002 
which is definitely larger than the value predicted in Q 
Q, by some 50 standard deviations. 

In an alternative scenario we could assume that the 
deviations from Eq. H10|l seen in Fig. ^Jd do not vanish in 
the limit — > cxd. In that case one should also allow for 
a modified scaling variable 



6 - (2/ - S)/{D + 1-25) 



(12) 



with an unknown small (non-universal) parameter 5. Val- 
ues of (£> -I- 1 - 26)p{y)/fp{^s) for D = 128 and three dif- 
ferent values of 6 are plotted in Fig. Similar results 
are obtained for other D. Combining them, we see that 
best scaling (i.e. least dependence of D for small y) is 
obtained for S « 0.02. For this value of S we have 

lim D'^/^p{y)/y^/^ = (0.95 ± 0.02) x A . (13) 

y — ^0, D — *oo 

The large uncertainty in this estimate reflects the rather 
steep slope of the central curve in Fig. ^ as ^5 — > and 
the related uncertainty in the best estimate of 5 If 
we accept Eq. (^21, we obtain B = 2.04 ± 0.04, which 
is in perfect agreement with 0, Q . Thus we have two 
scenarios. Both imply very large corrections to scaling. 
While the a priori simpler scenario would be in conflict 
with the theoretical prediction, this prediction suggests 
that indeed the second scenario is correct, which is our 
preferred solution. 

In Fig.|31 we plot the end-to-end distance per monomer 
{x) /N versus N for various widths D. These curves be- 
come horizontal as A'^ — *■ 00, i.e. (a:;) increases indeed 
linearly with N, limAr_,oo (x)/A^ = Ax. In order to find 
how Aa; scales with D, we plot it in Fig. El on a doubly 
logarithmic scale. As indicated by the dashed line, it is 
fitted perfectly by the theoretical prediction 



{x)/N ^ D^-^^" = D-^^^ . 



(14) 



We also estimated the density of wall contacts pb (num- 
ber of monomers at ?/ = 1 or at y = divided by 2{x)). 
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For each fixed value of D, this becomes independent of N 
as iV — > CX3. The asymptotic values, obtained from plots 
similar to Fig. are also plotted in Fig. El The full line 
corresponds to the very simple prediction [5j 

Pb - D-^ (15) 

which is independent of v and indeed holds also for Gaus- 
sian chains. Eq. (|15(l can be easily understood, in terms 
of the pressure exerted onto the wall: 

p,^p^^^ 15-1-1/^+1/3 = ^-2. (16) 

Finally, we show in Fig. [3 the distribution pend{y) of 
chain ends. We see that Pend{v) scales for large D (and 
for N ^ Z?!/", of course). Theoretically it is predicted 
that i23|| 

Pend{y)-V^''"' (17) 

for y <^ D. We see that this is indeed verified (the full 
lines in Fig. [3 correspond to the prediction). But in con- 
trast to the monomer density which was described fairly 
well for all y by the product of the power laws holding 
near the two walls, the same definitely does not hold for 
Pendiy)- There the function 5^(0 = const (^(1 -0)'^/^^ 
does not describe the behaviour away from the walls. 

IV. SUMMARY 

We have shown that we could simulate 2-d polymers, 
modelled as self-avoiding walks, with chain length up to 



125000 on strips of widths up to 320. This was possible 
using the PERM algorithm with Markovian anticipation. 
The fact that PERM gives by default very precise esti- 
mates of free energies allowed us to measure precisely 
the forces exerted onto the walls, by measuring how the 
critical fugacities depend on the width of the strips. We 
verified all critical scaling laws predicted for this problem, 
including the scaling of monomer and end point densities 
near the walls and the scaling of the total pressure with 
chain length and with strip width. 

The only prediction for which we found possibly dis- 
agreement is for the universal amplitude ratio B defined 
in Eq. (Q. A scenario based on some minimal assump- 
tion about scaling functions and corrections to scaling 
gives an estimate higher than the prediction by some fifty 
standard deviations. But a different scenario, maybe less 
plausible a priori but not very unlikely either, gives per- 
fect agreement with the prediction. This illustrates again 
that one should be very careful about corrections to scal- 
ing, and that even very precise simulations do not always 
give unique answers when their analysis is not guided by 
a reliable theory. 

Previous simulations of 3-d polymers between two par- 
allel planar walls had indicated that also there the value 
of B might be larger than predicted, but those simula- 
tions had very large uncertainties. Using PERM we can 
simulate fairly easily much longer chains with rather high 
statistics. Results of such 3-d simulations will be given 
elsewhere. 
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